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ABSTRACT 

Recent low-redshift observations have attempted to determine the star formation his- 
tories of ehiptical galaxies by tracking correlations between the stellar population 
parameters (age and metallicity) and the structural parameters that enter the funda- 
mental plane (size rg and velocity dispersion a). These studies have found that velocity 
dispersion, rather than effective radius or dynamical mass, is the main predictor of 
a galaxy's stellar age and metallicity. In this work, we apply an analytic model that 
predicts the structural properties of remnants formed in major mergers to progenitor 
disk galaxies with properties taken from two different semi-analytic models (SAMs). 
We predict the effective radius, velocity dispersion, luminosity, age, and metallicity of 
the merger remnants, enabling us to compare directly to observations of early- type 
galaxies. We note that our SAM population consists solely of elliptical galaxies whose 
evolution is artificially truncated following the formative major merger, while we are 
comparing to a low-redshift population of elliptical and SO galaxies that may have 
more complex formation histories. While we find a tight correlation between age and 
velocity dispersion, we find a stronger dependence of age and metallicity on effective 
radius than observations report. We find that these correlations arise as a result of the 
dependence of gas fraction, age, and metallicity on the stellar mass in the progenitor 
disk galaxies. These dependences induce a net rotation in the radius- velocity plane 
between the correlations with effective radius and circular velocity in the disk galaxy 
progenitors, and the correlations with effective radius and velocity dispersion in the 
elliptical galaxy remnants. The differences between our results and those from observa- 
tions suggest that major mergers alone are not sufficient to produce the observed lack 
of correlation between effective radius and stellar population parameters. Simulations 
have suggested that subsequent minor mergers may introduce scatter in the effective 
radius while leaving the velocity dispersion essentially unchanged. Incorporating such 
minor mergers into the model may, then, bring the simulations into closer agreement 
with observations. 

Key words: galaxies: interactions - galaxies: evolution - galaxies: elliptical and 
lenticular, cD - galaxies: formation 
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Simulations of galaxy formation have shown that ellipti- 
cal galaxies may be created through major mergers of disk 
galaxies (e.g. [Toomre k Toomre||1972| |Toomre|p977| |Mi- 



hos Hernquist||1994 ). Small, compact ellipticals can form 



© 0000 RAS 



2 Porter et al. 



at high redshifts from gas-rich 'wet' mergers (Dekel & Cox 



2006 Ciotti et al. 2007). Subsequent dissipationless 'dry' 
mergers greatly increase the galaxies' stellar masses and 
radii, producing the massive ellipticals seen in the local uni- 



verse (Khochfar & Burkert 


2003| |Naab et al. 2006 Naab 


leTal. 2009; Oser et al.|2012 


). 



Early type galaxies can be described by a two- 
dimensional plane relating effective radius (i?e), central stel- 
lar velocity dispersion (cr), and effective surface brightness 
(/e), termed the fundamental pl ane (FP) (|Djo rgovski fc| 



on effective radius. G09 posited that the lack of dependence 
on Re occurs because effective radius is strongly dependent 
on the orbital parameters of the major merger and subse- 
quent dissipationless mergers. This dependence introduces a 
large amount of scatter in the remnants' radii, effectively di- 
luting any corresponding radial relations that were present 
in the progenitors. 

Semi-analytic models (SAMs) provide a framework to 
simulate the formation and evolution of elliptical galaxies 



more, the fundamental plane is not an exact relation; galax- 
ies have a degree of scatter around the FP, in effect making 
the fundamental plane 'thick'. Observations indicate that 
this scatter increases with redshift, particularly among less 
massive galaxies ( Treu et ah]! 2005b ) . 



Davis||1987 Dressier et al. 1987 Faber et al. 1987). This 


et al. 19941 |Somerville & Primack |1999| |Cole et al.||2000| 


plane is tilted from the plane one would expect from a sim- 
ple application of the virial theorem, indicating that further 
processes, such as non-homology or a varying mass-to-light 


Hatton et al. |2003| |Croton et al. |2006 


De Lucia et al. 


20061 IBower et al."2006; 'Somerville et al.l 


2008b^ "Fontanot 


et al.t2009i Benson & Bower 2010; Cook et al. 2010; Guo 


ratio, must have an effect (|j0rgensen et al.||l996|). Further- 


et al.||2011| |Somerville et al.||201l|. However, attempts to 



More specifically, while the slope of the FP appears 
unchanged for high-mass ellipticals since z ^ 1, low-mass 
ellipticals at high redshifts have higher surface brightnesses 
than their effective radii and velocity dispersions would seem 



to predict ( 


van der Wei et al. 2004| Treu et al. 


2005a|b 


J0rgensen et al.||2006 van Dokkum & van der Marel||2007). 



If we consider a projected FP, where surface brightness is the 
dependent parameter, then these low-mass galaxies tend to 
lie above the mean FP relation at high redshift, and to fall 
onto the FP over time. 

There are indications that this residual thickness in the 
FP correlates with the stellar population age. Forbes et al.| 
( 1998| ) and [Terlevich k Forbes] ( |2002| ) found that galaxies 



with higher residual surface brightnesses are younger than 
those that lie near the mid-plane of the FP; conversely, those 
with lower residual surface brightnesses are older. 

More recently, observations from the Sloan Digital Sky 
Survey (SDSS) have been used to analyze stellar popu- 
lation trends both within the R-cr projection of the FP, 
and through the thickness of the FP, using residual surface 
brightnesses ("Gr aves et aLl|2009a|b[ |2010 Graves & Faber 
j2010). By stacking spectra of galaxies with similar stellar 
properties and measuring the Lick indices on those spec- 
tra, the authors were able to derive [Fe/H],[Mg/H], [Mg/Fe], 
and stellar age for a population of passive early-type galax- 
ies. In agreement with [Forbes et al.| ( |19 98) and Terle vich k\ 
[Forbes (2002), Graves et al. (2009b, hereafter G09) found 
that younger galaxies lie above the FP, having relatively 
higher surface brightnesses, while older galaxies lie below it; 
galaxies above the FP also tended to have higher [Fe/H] and 
[Mg/H], and lower [Mg/Fe]. G09 also determined that age, 
[Fe/H], [Mg/H], and [Mg/Fe] increase with velocity disper- 
sion throughout the FP, independent of the residual surface 
brightness. These same properties are almost independent of 
i?e, indicating that a galaxy's velocity dispersion is a more 
reliable determinant of its star formation history than its dy- 
namical mass (oc a^Re). The strong dependence of age and 
metallicity on velocity dispersion is consistent with previous 
studies (Smi th et al.|2007|[Neian et aLpOOSf . Similar results 
were obtained in a recent analysis of the Six-degree Field 



Galaxy Survey (Jones et al.||2004| 2009) by Springob et al. 



use SAMs to study the fundamental plane have been lim- 
ited by the difficulty in modeling the effective radii and ve- 
locity dispersions of ellipticals. Cole et al. (2000) provide 



a simple formula to predict the radii of elliptical galaxies 
following a major merger using the virial theorem and con- 
servation of energy. While this relation may be correct for 
dissipationless gas-poor mergers, the energy lost due to star 
formation in gas-rich mergers results in a deviation from 
the virial relation, and smaller remnant radii. Incorporating 
this dissipation is probably essential: a recent study using 



the Bower et al. (2006) SAM framework could match the 
observed SDSS size-mass relation only by adding dissipa- 
tion (Shankar et al.|2011 ). A second recent SAM paper that 
did not include dissipation in modeling merger remnants had 
a too shallow slope (Guo et al. 201 1|). The predicted scat- 
ter in the size-mass relations in these papers was also much 
larger than observed. Full N-body/SPH simulations of gas- 
rich mergers have bee n more successful ([D ekel & C ox|2006| 
Robertson et al.|2006 Hopkins et al.|20Q8| , but they do not 
simulate enough galaxies to compare with observations of 
populations of galaxies as in G09. 

Using hydrodynamical simulations ( Cox 2004] jCox et al 



2006| |2008t , Covington et al. (2008, hereafter C08) devel- 
oped an analytic model to predict the effective radius and 
velocity dispersion following the major merger of two disk 
galaxies. This model was based upon the virial theorem and 
additionally incorporated energy losses due to dissipation. In 
further work, Covington et al. (2011, hereafter Cll) simpli- 
fied the model and applied it to mergers of disk-dominated 
galaxies with properties taken from the Croton e t al.[ ( [2006| 



Millennium SAM and the Somervill e et al.[ ([20081o) SAM 



via post-processing. Cll showed that the C08 model repro- 
duces the observed steepening in the size-mass relation of 
ellipticals when compared to disks, as well as the observed 
decreased dispersion in radius compared to the progenitor 
galaxies (Shen et al.|[2003| . In addition, this model quali- 
tatively reproduced the evolution of these properties versus 



redshift (Trujillo et al. 2006 ). The C08 model was also shown 
to correctly reproduce a tilt in the FP away from the simple 
virial relation. Using the methods of C08 and Cll, as well 
as an alternative prescription from [Hopkins et al. (2009), 
Shankar et al. (2011) reached similar conclusions. 



(2011), though this work finds a slightly higher dependence 



The predictions of effective radius and velocity disper- 
sion from C08 and Cll, combined with information from 
SAMs, allow us to determine a simulated galaxy's location 
above, within, or below the FP, as measured by its resid- 
ual surface brightness. In this work, we relate the age and 
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metallicity of simulated galaxies across the FP with effective 
radius and velocity dispersion. We compare galaxies with 



initial conditions taken from the Croton et al. (2006) and 



Somerville et al. ( 2008b ) SAMs to each other as well as to the 



analysis of G09, and find that the stellar population metalli- 
cally is correlated with both radius and velocity dispersion, 
while age is strongly dependent on velocity dispersion and 
almost independent of radius. The metallicity correlation is 
in contrast to observational results (G09), which found a 
tight correlation with velocity dispersion only and almost 
no correlation with radius for either age or metallicity. We 
find that all five measured parameters (age, metallicity, lu- 
minosity, velocity dispersion and radius) are closely tied to 
those of the disk progenitors, with the main difference be- 
ing a rotation in the two-dimensional Re-cr projection of the 
fundamental plane due to differing amounts of dissipation 
in galaxies with different gas fractions. We find no change 
in the FP-stellar population correlations with the redshift of 
the formative merger. Examining the thickness of the fun- 
damental plane, we find that galaxies that fall below the FP, 
with relatively lower surface brightnesses, tend to be both 
older and more metal-poor than galaxies that lie above the 



FP. This is in agreement with observations ([Graves et al 



2009b Terlevich & Forbes 2002 i Forbes et al. 1998) 



Section [2?T] describes the two semi- analytic models used 
in our analysis. Section [2.2| provides an overview of the an- 
alytic model used to calculate the radius and velocity dis- 
persion of the remnants of major mergers between two disk 
galaxies. The prescription we use combines the more de- 
tailed model for effective radius and velocity dispersion from 
C08 with an improved model to predict the concentration of 
dark matter within the effective radius Re ( [Covington et al.| 
2011). Sections 2.3 through 2.4 describe our calculations of 



the light-weighted ages, light- weighted metallicities, and lu- 
minosities of the remnants at redshift zero, using the mass- 
weighted ages and metallicities of the progenitors, as well as 



the stellar population models of Bruzual & Chariot ( 2003 ) 



We also describe the procedure for separating galaxies ac- 
cording to their location on the fundamental plane, based 
on their residual surface brightnesses. Section 3 presents a 
summary of the low-redshift observations of G09, to which 
we make direct comparisons. 

We present results beginning in Section 4, in which we 
examine the relationships between age and metallicity, on 
the one hand, and velocity dispersion and radius, on the 
other, throughout the fundamental plane for the |Croton[ 
et al.| ([20061 and [Somerville et al.[ ([2008b|) SAMs. In both 



cases (age and metallicity), and in both SAMs, the relations 
for the remnants represent a net < 20° counterclockwise ro- 
tation from relations for the progenitors; thus, correlations 
with radius (velocity dispersion) in the progenitors tend to 
produce mixed correlations involving both components in 
the remnants. Section 5 examines the residual thickness of 
the FP and the structural differences between galaxies ly- 
ing above and below the FP. In Section 6 we compare our 
results to those of G09. Section 7 examines the evolution 
of the FP correlations as a function of the redshift of the 
formative merger event, and attempts to predict the effects 
of including subsequent evolution via multiple mergers, in- 
cluding minor mergers. 



2 METHODS 

We begin with an overview of the main components of the 



model: the Millennium ( [Croton et al.[[2006| ) and [Somerville 



et al. (2008b), hereafter S08, semi-analytic models, from 



which we drew progenitor galaxies, and the dissipational 
merger model, which was used to predict the structural 
properties of the merger remnants. 



2.1 The semi-analytic models 

We model elliptical galaxy formation using the properties of 
late- type progenitors, taken from SAMs, to predict the ages, 
metallicities, luminosities, velocity dispersions, and effective 
radii of elliptical remnants following a major merger. Us- 
ing the results of N-body dark matter simulations or the 
extended Press- Schechter formalism, the two SAMs con- 
struct merger trees of dark matter halos, which are then 
seeded with galaxies. When two halos merge, their galax- 
ies are considered to merge on a characteristic timescale, so 
long as the smaller galaxy is not tidally disrupted before 
the merger completes. The evolution of this galactic sub- 
structure is treated in a semi-analytic manner, using sim- 
ple prescriptions for photoionization, radiative cooling, star 
formation, supernova feedback, and black hole growth and 
feedback. This methodology has been shown to accurately 
reproduce low-redshift observed relations, such as the stellar 
mass function, cold gas fractions, and specific star formation 
rates. For more detail, we refer the reader to [Somerville et ah] 
( [2008b| ; [Croton et al.[ ((2006). 

Drawing from the Millennium and S08 merger trees, we 
select pairs of disk-dominated galaxies that merge with a 
stellar mass ratio equal to or greater than 1:3 and with a 
stellar mass greater than 10^ M0. We exclude mergers of 



bulge-dominated galaxies, since the SAMs do not currently 
produce accurate estimations of bulge size, and hence stellar 
radius, for these galaxies. These mergers also tend to be dis- 
sipationless, as the bulge-dominated galaxy is usually both 
more massive and gas-poor; thus the baryonic gas fraction of 
the merger is relatively low. Since we want to produce quies- 
cent galaxies at redshift zero, we exclude any mergers from 
the past 1 Gyr (z < 0.075). Making these selection cuts, we 
include 57,435 mergers from the S08 SAM and 53,416 merg- 
ers from the Millennium SAM. For simplicity, we do not 
consider any further evolution following the major merger; 
possible implications of subsequent evolution are discussed 
in section 7.1. We also do not attempt to model the for- 
mation and growth of bulges via disk instabilities, or the 
formation of SO galaxies. In effect, we are only studying the 
relations elliptical galaxies would have if they formed in a 
major merger and then passively evolved to redshift zero. 

Both simulations use cosmological parameters taken 
from the first-year Wilkinson Microwave Anisotropy Probe 
(WM AP) results: Qm = 0.25, = 0.045, h = 0.73, ag = 
0.90 ( jSpergel et aL][2U03 ). 



2.2 Dissipational model for remnant structural 
parameters 

To locate elliptical remnants in FP-space, we have built 



upon the dissipational model of Covington et al. ( Covington 
et al.|2008 Covington et al.|2011 ). This is a simple analytic 
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model to predict the effective radius and velocity dispersion 
following major mergers, derived from adding radiative en- 



ergy loss to a standard dissipationless model (see Cole et al 



2000). Here we provide a brief summary of the derivations 



used to calculate the parameters we will be using, namely 
the remnant galaxy's stellar mass, half-mass radius, and ve- 
locity dispersion. For a full description, we refer the reader 
to COS and Cll. 

In a major merger, the efficiency of the burst of star 
formation depends on the amount of cold gas available and 
the dissipation due to the merger. We parameterize the dis- 
sipational strength of each galaxy as fk = AE/Ktot, where 
AE is the impulse between the two galaxies and Ktot is 
the total initial kinetic energy of the galaxy. To determine 
the impulse, pericentric merger distances and velocities were 
drawn from the orbital parameter distributions given in |Ben-| 
( 2005 ) , though the results do not change significantly if 



newer parameters fWetzel |2010 ) are used. The mass of new 
stars formed for each progenitor is given by: 



-^^new stars — C new gas fh 



(1) 



where Mgas is the mass of gas in the progenitor and Cnew 
0.3 is obtained by fitting to the "Cox" (2004 ) simulations. 

Equating the initial internal energies of the progenitors 
with the sum of the internal energy of the remnant and the 
radiated energy lost due to star formation, COS found: 



-E-init — CintG 



(Ms,l +Mns,l)^ , (Ms,2 + Mns,2 



Ri 



R2 



Efinal — ^init + ^rad , 

and 

_^ ^/(Ms,l +Ms,2 + Mns,l +Mns,2)^ 
^final — OintCrl 



i^final 



(2) 



(3) 



(4) 



where Ms,i is the stellar mass of the progenitor, Mns,i is the 
mass of stars formed in the progenitor during the merger, Ri 
is the stellar half-mass radius of the progenitor, and dnt = 
0.5 is a constant tuned to the Cox| ( |2004J simulations. 
The radiated energy is given by: 



^rad = iri/g,i/k,i(l + fk,i) 



(5) 



where Ki, /g,i, and /k,i are the total kinetic energy, baryonic 
gas fraction, and fractional impulse of progenitor i. An addi- 
tional term, corresponding to the energy of the orbit itself, 
is added to the conservation of energy calculation (equa- 
tion 2.2); however this term is typically on the order of a 



few percent of the radiated and initial energy terms, and 
so does not significantly affect the remnant's radius. Using 
equations |2.2| and |2.2| the half-mass radius of the remnant 
can be determined. 

We use a modified virial-type relation to determine the 
velocity dispersion of the remnant: 

where Ms,f is the stellar mass of the remnant, /dm,f is the 
central dark matter fraction of the remnant (i.e. the fraction 
of dark matter within the stellar half-mass radius), Ri is the 
half- mass radius of the remnant, and Csie- = 0.15 is a con- 



energy terms of equations |2 . 2| [2^2] and [2T2] use the total mass 
of the progenitors, rather than the baryonic mass. 

2.3 Calculating the light-weighted remnant age 
and metallicity 

As we interested in the effects of major mergers alone, we 
artificially freeze all the physical properties of the galaxy di- 
rectly after the merger. The stellar population is allowed to 
passively evolve, but we do not consider subsequent mergers, 
star formation, or mass accretion. Both the SOS and Millen- 
nium SAMs provide mass-weighted ages and total metallic- 
ities of the progenitor galaxies at the time of the merger. 
In order to directly compare to the G09 data, we convert 
the mass- weighted values to light- weighted ones, consider- 
ing the remnant galaxy to be comprised of two progenitors 
and a new stellar population formed in a burst during the 
merger. The weighting formula we have used is: 

S(Am,p * Lp) + An * ivN 



(7) 



S(Lp) 

where Al,r is the light- weighted age of the remnant, Am,p 
is the mass- weighted age of the progenitor, An is the age 
of the new stars, and Lp and Ln are the V-band luminosi- 
ties of the progenitor and new stars, respectively. A similar 
formula is used to determine the light-weighted metallicity, 
where the metallicities of the progenitors and new stars are 
used in place of their ages. The summations are over pro- 
genitor galaxies. All ages and luminosities are calculated at 
a redshift of 0.0, to compare with G09. 

Since the stellar populations of the progenitor galaxies 
are assumed to passively evolve to redshift zero, we take 
their redshift zero mass-weighted ages to be equivalent to 
their light-weighted ages, setting the mass-to-light ratio to 
1. We assume all new stars form at the instant of the merger, 
so the age of the new stars at redshift zero is equivalent to 
the lookback time of the merger. 



Following [Croton et al.| ( [2006] ) and Somerville et al 
(|200Sb'), we assume that stars formed during the merger 



have the same metallacity as the cold gas of their host 
galaxies at the onset of the merger. While the semi-analytic 
models include prescriptions to model chemical evolution 
through supernova feedback and accretion from the inter- 
galactic medium, we do not include any further chemical 
evolution after the merger event. 

2.4 Computing the remnant luminosities 

The luminosities of the progenitors and remnants are de- 
termined using the stellar population synthesis models of 



[Bruzual Chariot ( 2003| , hereafter BC03. As noted above, 
all calculations are done at a redshift of zero. 

We use a model based on the Padova 1994 stellar evo- 



lution tracks (Bertelh et al. 1994) with a Chabrier (2003) 



stant set to match the Cox (2004) simulations. The internal 



initial mass function. The progenitor galaxies are taken to 
be simple stellar populations (SSPs), with luminosities at 
the time of the merger provided from the respective SAMs. 
We use these quantities to calculate the mass-to-light ratio, 
and hence the luminosity, at redshift zero using a linear in- 
terpolation of the BC03 models in log time and (roughly) 
log metallicity. We exclude progenitors whose metallicities 
are outside the range of the BC03 models, 0.0001 ^ Z ^ 
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0.05; these galaxies represent less than 0.3% of the progeni- 
tor population and would not significantly affect the results. 

The calculation of the mass-to-light ratio for new stars 
formed during the merger is slightly more involved. The star 
formation rate from the merger is modeled as exponentially 
decaying in time, and we use BC03 to form tracks based on 
the Padova 1994 evolution tracks for this population. More 
specifically, the star formation rate is: 

[IMo - /rMpgWJt-' exp(-t/T) 



(8) 



where MpG(t) is the mass of processed gas that has been 
converted into stars and returned to the interstellar medium, 
determined within BC03, and /r is the fraction of Mpg 
recycled back into new star formation, r is the e- folding 
timescale; for the actual calculation we used an equiva- 
lent parameter, the efficiency of star formation: e = 1 — 
exp( — 1 Gyr/r). Note that BC03 defines efficiency as the 
fraction of gas converted into stars after 1 Gyr, while the 
efficiency we use is the fraction of gas converted over the 
entire burst. We use /r = 0.43, which is appropriate for the 
adopted Chabrier IMF. 

2.5 Binning in the fundamental plane 

Once the age, metallicity, effective radius, luminosity, and 
velocity dispersion have been calculated for the remnants, 
we separate them into five slices according to their loca- 
tion above or below the fundamental plane, using surface 
brightness as the independent variable. Since we intend to 
compare to the G09 results, only galaxies that fall within 
the G09 range of radius, (0.0 kpc < log Re < 0.7 kpc), and 
velocity dispersion, (2.0 km s~^ < log a < 2.4 km s~^), are 
included in the fitting routine. We use a least-squares fit to 
determine a relation between (log) Re, (log) a, and (log) /e, 
finding 



log le = -1.28 log Re + 1.32 log a - 0.29 
for SOS and 



log le 



-1.31 log Re + 1.22 log a - 0.29 



(9) 



(10) 



for Millennium. For each remnant, the predicted sur- 
face brightness is determined using the above rela- 
tion, and galaxies are separated by their residuals, 
where (-0.3 L© pc"^ < A log le < -0.1 L© pc"^), 

(-0.1 L© pc~^ < A log le < 0.1 L© pc~^), and 
0.1 L© pc"^ < A log le < 0.3 L© pc ^) are 
termed the low-, mid-, and high-FP, respec- 
tively. Galaxies with residuals outside the range 
(-0.3 L© pc"^ < A log le < 0.3 L© pc ^) are placed 
in the very low ('vlow') and very high ('vhigh') FP slices. 

If we plot the galaxies according to their location in 
FP-space, 97% (93%) of the SOS (Millennium) galaxies fall 
within the low-to-high FP slices, with the majority of the 
outliers falling below the FP (Figure [T]). 

After separating the remnants by their location within 
the FP, we place them in bins according to their radius and 
velocity dispersion. We then calculate the median age and 
metallicity for all galaxies within each bin, excluding bins 
with fewer than five galaxies for statistical purposes. These 
values are used to form contours relating the stellar pop- 
ulation parameters, namely age and metallicity, with the 
fundamental plane parameters and residuals. 



Comparing the SOS and Millennium populations, we 
find that they occupy slightly different regions of the Re-cr 
parameter space (Figure [2]). The Millennium galaxies tend 
to have larger radii and smaller velocity dispersions than 
the SOS galaxies. This is a result of the differing progeni- 
tor populations, as noted in Cll; Millennium disk galaxies 
are larger and more gas-poor than SOS disk galaxies. As 
a result, major mergers induce less dissipation, producing 
remnants with larger radii. Since the velocity dispersion is 
inversely proportional to the radius in the model, the Mil- 
lennium galaxies have smaller velocity dispersions as well. 
As a reference, we have overplotted the region analyzed in 
G09, seen as the black boxes in the middle three FP slices, 
as well as contours demonstrating the G09 parameter space. 
Both SAMs display a population of galaxies with low radii 
that is not seen in G09. These galaxies tend to have the the 
highest merger redshifts, suggesting that they correspond 
to compact ellipticals at high redshifts and are not likely to 
be seen in large numbers in the local universe. The Millen- 
nium SAM also contains a population of galaxies with low 
velocity dispersions, log a < l.S km s~^, which falls below 
the completeness threshold of G09. It is worth noting that 
the Millennium parameter space roughly resembles that of 
G09, with the aforementioned higher radii and lower velocity 
dispersions. 



3 SUMMARY OF OBSERVATIONS 

We compare our findings to a recent survey of early-type 



galaxies from the Sloan Digital Sky Survey (SDSS) (York 
2000) Spectroscopic Ma in Galaxy Survey ([Strauss 



et al. 



et al. 



2002| ) Data Release 6 ( [Adelman- McCarthy et al 



200S). The sample of galaxies is described in Graves et al. 
( |2009a|bt . To summarize, the galaxies selected were ob- 
served in the redshift range 0.04 < z < O.OS, with light 
profiles that were both centrally concentrated and fit a de 
Vaucouleurs profile. To prevent a small proportion of young 
stars from biasing the measured luminosity, G09 excluded 
actively star- forming galaxies. Using colors and emission-line 
intensities, G09 also rejected Seyfert hosts, low ionization 
nuclear emission-line region (LINER) hosts, and transition 
objects, as they can host active galactic nuclei (AGN) which 
have been found to have light profiles intermediate between 
early- and late- types ( [Kauffmann et al.|2003| ) . Although we 
have excluded actively star-forming galaxies from the sam- 
ple, as noted above, no attempt has been made to exclude 
galaxies that the SAMs characterize as having active AGN 
at redshift zero. 

We also note that the population of simulated galax- 
ies consists solely of ellipticals that have formed in a ma- 
jor merger with no subsequent evolution, while the G09 
population may contain a significant number of passive SO 
and Sa galaxies, as well as ellipticals that may have formed 
through alternative processes or undergone structural evolu- 
tion following their formation. An analysis of quiescent red 
sequence SDSS galaxies similar to the population studied 
in G09 found that 36% of the galaxies were true elliptical 
galaxies, 15% were SO galaxies, and 4S% were Sa galaxies 
(Cheng et al.||2011| ). Furthermore, while both the SOS and 
Millennium SAMs reproduce the early-type mass function 
when galaxies are separated by their stellar bulge-to-total 
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Figure 1. Distribution of galaxies through the fundamental plane for SOS (left) and Millennium (right). Galaxies are fit to a linear 
relation (horizontal axis) relating surface brightness with velocity dispersion and radius. The measured surface brightnesses are then 
plotted against the expected values. The areas between the solid black lines represent the slices we term the 'low-FP', 'mid-FP', and 
'high-FP', from bottom to top, according to the residual in surface brightness. Each slice has a thickness of 0.2. The areas below and 
above the red lines represent the 'vlow-FP' and 'vhigh-FP', respectively. 97% and 92% of the SOS and Millennium galaxies, respectively, 
fall within the middle three FP slices. The red, green, and blue contours enclose 50%, 70%, and 90% of all galaxies, while the grey points 
represent individual galaxies. 
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Figure 2. Distribution of radius and velocity dispersion for galaxies within each slice of the FP for SOS (top). Millennium (middle), 
and G09 (bottom). From left to right, the panels represent the 'vlow-FP', 'low-FP', 'mid-FP', 'high-FP', and 'vhigh-FP' slices. The grid 
lines show the bin definitions in the region observed in G09; the median age and metallicities are calculated within each bin. Compared 
to the G09 population, the Millennium galaxies are more likely to have low velocity dispersions, owing to the effects of incompleteness 
in the observations. Both SAMs contain a population of galaxies with low radii and high merger redshifts; these correspond to massive, 
compact ellipticals and are not likely to be seen in the local universe. The red, green, and blue contours enclose 50%, 70%, and 90% of 
all galaxies, while the grey points represent individual galaxies. 
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ratios, the selection cuts we have made in this paper hmit 
the sample to ^ 1/3 of bulge-dominated galaxies. Thus, a 
detailed comparison between the simulated galaxies and ob- 
servations will require either limiting the G09 sample to pure 
elliptical galaxies, or expanding the simulation to include a 
broader range of red sequence galaxies. This analysis will be 
performed in future work (Porter et al., in prep.). 

Using the Lick indices f Wort hey et al. 1994 Wort hey 
& Ottaviani |1997| on 16,000 stacked spectra, G09 cal- 
culated mean light-weighted ages and metallicities in 
bins with residual surface brightness above, within, and 
below the fundamental plane. The bins covered the 
approximate range (1.9 km s~^ < log a < 2.4 km s~^), 
(0.0 kpc < log Re < 0.7 kpc), and 
(-0.3 L© pc~^ < A log le < 0.3 Lq pc~^), where A log h 
is the residual surface brightness resulting from a log fit in 
radius and velocity dispersion. 



log le = 1.16 log a - 1.21 log Re + 0.55. 



(11) 



G09 formed contours relating the mean light-weighted 
age and light- weighted metallicities, [Fe/H], [Mg/H], and 
[Mg/Fe], to effective radius and velocity dispersion across 
three slices of the fundamental plane. The Millennium and 
SOS SAMs measure the total metallicity Z but do not sepa- 
rate contributions from different elements and only include 
Type II supernova; thus, we consider the total metallicity 
from SAMs to be most similar to [Mg/H], a measure of 
alpha-enhancement. These results can be seen in Figures 
7 and 9 of G09. The authors found that stellar population 
age and metallicity are nearly independent of effective ra- 
dius but strongly correlated with velocity dispersion. In all 
three slices of the FP, galaxies with larger a had older ages 
and higher metallicities. Stellar population age was also in- 
versely correlated with residual surface brightness, so that 
the youngest galaxies tend to fall above the FP, in agree- 
ment with earlier observations ( Forbes et al.|1998 Terlevich 
Forbes|[2Q02| ). 



4 RESULTS 

Our major findings are the contours seen in Figures [3] and 
|4] For both SOS and Millennium, stellar population age in- 
creases strongly with velocity dispersion and decreases with 
radius. We note that our parameter space has a much larger 
range in radius than velocity dispersion, so that while the 
contours appear nearly vertical, the radial dependence is 
non- negligible. Looking across the FP, galaxies that lie above 
the FP (those with the largest residuals in log le) have 



younger ages, as suggested in Graves et al. 



etaL] ( ri99SJ ; and [Terlevich Forbes] ( |2002t . As the Millen- 



(2009b); Forbes 



nium galaxies sample a larger portion of the 2-dimensional 
Re — (J parameter space (see Figure [2| , galaxies from the 
Millennium SAM tend to exhibit a larger gradient in age 
within a slice of the FP. 

If we compare the metallicity contours, the differences 
between the SAMs are slightly more pronounced. In both 
cases in Figure 4, metallicity increases with velocity disper- 
sion and radius. However, the Millennium remnants have 
a stronger correlation with radius, while for SOS the trends 
with radius and velocity dispersion are roughly equal. Galax- 



ies that lie above the FP tend to have higher metallicities, 
as in G09. 

To better quantify the dependence of the stellar pop- 
ulation parameters on the FP parameters, we have plotted 
the binned stellar age and metallicity as a function of ra- 
dius and velocity dispersion for the midplane of both SAMs 
(Figure |5|. We have also calculated the least-squares fit and 
Spearman rank coefficients for these relations. These two- 
dimensional relations confirm the trends seen in the con- 
tours. Both SAMs have a stronger dependence on velocity 
dispersion (p = 0.74 and 0.72 for SOS and Millennium, re- 
spectively) than on effective radius (p = -0.56 and -0.45 
for SOS and Millennium, respectively) with regards to stel- 
lar age. The Millennium metallicities have a nearly perfect 
correlation with radius (p = 0.96) and a weak correlation 
with velocity dispersion (p = 0.39), while the SOS metal- 
licities have a somewhat stronger correlation with radius 
(p = O.Sl) than with velocity dispersion (p = 0.60). With 
the exception of the Millennium metallicity-radius relation, 
all of the correlations have a large degree of scatter around 
the least-squares fit. 

There are marked differences between the ranges of 
metallicities in SOS and Millennium. Comparing the 25th 
and 75th percentiles, the range of metallicites in SOS 
(-0.50 < logZ < -0.24) is roughly twice that of Millen- 
nium (—0.26 < logZ < —0.14). While the average metal- 
licity in SOS is lower than in Millennium (logZ = —0.37 
and —0.20, respectively), the overall gradient in metallicity 
is also much larger in SOS, with the low (high) a galaxies 
having lower (higher) metallicities than in Millennium. 

We defer a detailed comparison between the two SAMs 
to future work, but we note here that the two SAMs use 
similar prescriptions for metal enrichment, modeled after 



De Lucia et al. (2004); for every solar mass of star forma- 



tion, a fixed yield of metals is returned to the cold gas disk. 
There are however key differences between the two SAMs. 
The SOS progenitor galaxies span a larger range in metallic- 
ity as a function of stellar mass (see Figure 11 right pane). 



Furthermore, the average stellar mass of a remnant is lower 
in SOS than Millennium, with values of log Mq — 10.37 and 
10.52, respectively. Finally, the SOS galaxies tend to merge 
at higher redshifts than their Millennium counterparts, with 
the average merger redshift being 1.73 and 0.99, respectively. 



2011 



Since metallicity tends to increase with stellar mass (Gal- 
lazzi et aL]|2005 ) and decrease with redshift (Pipino et al. 
""^^■'Laskar et ar20iri 'Calura et al. 2009] [Maiohno et~ 

6}, 



200S| iPipino et al. 200S, Erb et al. 2006t , the latter two 
trends bias the SOS galaxies toward lower metallicities. How- 
ever we note that, even at fixed redshift, the mass- metallicity 
has a larger gradient in SOS than in Millennium (see Figure 
right panel). 

If we neglect the range in metallicity within the SAMs, 
we find that the two SAMs produce remarkably similar age- 
FP and metallicity-FP relations. As we will explain in more 
detail later, this agreement stems from the correlations be- 
tween the progenitor gas fraction, age, and metallicity, on 
the one hand, and the progenitor stellar mass, on the other. 
Since the progenitor populations of the two SAMs have simi- 
lar stellar mass dependences, the remnants have similar fun- 
damental plane relations. 

Stellar population trends through the thickness the FP 
can arise in a number of different ways; any process that in- 
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Figure 3. Relation between luminosity- weighted age, effective radius, and velocity dispersion for elliptical galaxies in SOS (top) and 
Millennium (bottom). The different panels represent slices of the FP, as shown in Figure[2] In both SAMs, stellar population age increases 
strongly with velocity dispersion and decreases with radius. Galaxies that lie above the FP also tend to be younger than those that lie 
below the FP. These results are in qualitative agreement with observations, with the exception of the dependence on radius. 



creases the dynamical mass-to-light ratio of a galaxy will 
move it further below the virial plane. [Graves Faberl 
(2010) provide a decomposition of the deviation from the 



virial theorem, separating it into four components: 

(i) The ratio of the estimated dynamical mass to the true 
mass within one Re- 

(ii) The ratio of the true mass within one Re to the pro- 
jected stellar mass. 

(iii) The ratio of the projected stellar mass with an as- 
sumed initial mass function (IMF) to the true projected stel- 
lar mass and its corresponding IMF. 

(iv) The stellar mass-to- light ratio for the assumed IMF. 

For the simulated galaxies and their corresponding FP, the 
first and third terms are identically one, as we have no un- 
certainty in the dynamical mass estimate and we model and 
'observe' galaxies using the same IMF. Since we know the 
stellar mass-to-light ratio for the galaxies and can calculate 
the central dark matter fraction (DMF, see below), we may 
calculate the second and fourth terms directly. 

We have used the same process as described above to 
project the central dark matter fractions and stellar mass- 
to-light ratios of elliptical galaxies across and through the 
FP (Figures [g] an d [7| . The result s are in agreement with 
the conclusions of Graves & Faber ("2010 ): galaxies that fall 



radius, while the mass-to-light ratio increases with velocity 
dispersion. The former result is not surprising, as the ra- 
tio of dark matter to baryonic matter increases with radius; 
thus, larger galaxies enclose more dark matter. Comparing 
the trends through the thickness of the FP, at fixed Re and 
a the dark matter fraction has a larger degree of variation 
than the mass-to-light ratio. If we limit the sample to the 
middle three FP planes analyzed in G09 and [Graves Faber] 
( 2010), we find that the average variance in the dark matter 
fraction contributes 80% and 88% of the thickness of the 
FP for S08 and Millennium, respectively. This is in general 
agreement with Graves & Faber (2010), who found that the 



below (above) the FP have higher (lower) dark matter frac- 
tions and higher (lower) stellar mass-to-light ratios. Fur- 
thermore, the dark matter fraction tends to increase with 



dark matter fraction and IMF variation have a combined 
contribution in the range of 47% - 98%. This is an indi- 
cation that underlying structural differences, as opposed to 
stellar population effects, are the main contributors to the 
thickness of the FP. 

The variation in the central dark matter fraction can 
arise in a variety of ways: an increase in the stellar mass 
within one effective radius, an increase in the concentration 
of the host halo, and an increase in the effective radius at 
a given stellar mass will all have the effect of increasing the 
central dark matter fraction. In the dissipational model, the 
central dark matter fraction is estimated by calculating the 
halo half-mass radius, assuming that the two progenitor ha- 
los merge dissipationlessly, and fitting the resulting radius 
to an isothermal profile to calculate the mass of dark matter 
within one stellar effective radius. At fixed mass and stellar 
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Figure 4. Relation between luminosity-weighted metallicity, effective radius, and velocity dispersion for elliptical galaxies in SOS (top) 
and Millennium (bottom). The different panels represent slices of the FP, as shown in Figure [2] In both SAMs, metallicity increases with 
radius and velocity dispersion (contours are tilted); this is in contrast to observations, in which metallicity increases solely with velocity 
dispersion. Both SAMs display a weak trend through the FP, such that galaxies that lie above the FP have slightly higher metallicities, 
in agreement with observations. Galaxies from the Millennium SAM span a much lower range in metallicity. 



half-mass radius, halos with larger dark matter half-mass 
radii will have lower concentrations and lower central dark 
matter fractions. We choose an isothermal profile, as both 
simulations and observations have suggested that the central 
regions of dark matter halos containing early-type galaxies 
are isothermal in nature ( Auger et al.|2010b|a Bolton et al. 



2008| [Duffy et al.||2010| [Koopmans et al.||2066| ). For a dissi- 
pationless merger, 



m,f 



(Mdm,l + Mdm,2) 
+ 



(12) 



where Mdm,i and Mdm,2 are the total masses of the dark 
matter haloes of the two progenitors, Rdm,i and Rdm,2 are 
the half-mass radii of the haloes of the progenitors, and 
^dm,f is the halo half-mass radius of the remnant. Progeni- 
tors with small halo half-mass radii will produce a remnant 
with a correspondingly contracted halo and a higher concen- 
tration of dark matter within one Re. Thus, the variation of 
the dark matter fraction at fixed Re through the FP arises 
from variations in the concentrations of the progenitor spiral 
galaxies. 

The variation in the stellar mass-to-light ratio can fur- 
ther be separated into two components. Galaxies that merge 
at lower redshifts will have more recent star formation and 
decreased mass-to-light ratios, but there may also be an in- 
trinsic variation in galaxies at a fixed redshift due to differ- 
ences in stellar age. When we calculate the median mass-to- 



light ratio and merger redshift within each bin in Re and a, 
and compare these values through the FP, we find no cor- 
relation between merger redshift and stellar mass-to-light 
ratio in the lower four FP slices, while the vast majority of 
the galaxies in the highest FP slice have low merger redshifts 
and low mass-to- light-ratios (Figure [sj. This suggests that 
the effects of recent mergers are limited to galaxies that lie 
far above the FP, and that variations within the other four 
slices represent intrinsic differences. This is consistent with 



Graves & Faber (2010), who found that the thickness of the 



FP was due to structural differences in the stellar mass sur- 
face density rather than the passive fading of galaxies with 
young stars. 



4.1 Role of progenitor scaling relations 

Since our work is based on merging late-type galaxies with 
properties taken from the Millennium and SOS SAMs, the 
results will reflect any biases present in the SAM popula- 
tions. Comparing the evolution of the Millennium progen- 
itors' size-mass relations to observational data, Gil found 
that spiral galaxies that merged at the lowest redshifts, 0.0 
< z < 0.5, were on average 50% too large, while the most 
massive galaxies were too large at all redshifts. The slope 
of the late-type size-mass relation, however, was consistent 
with observational results, both in the local universe and at 
high redshifts. The size-mass relation for elliptical galaxies 
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Figure 5. Dependence of stellar age and metallicity on radius and velocity dispersion for elliptical galaxies in SOS (left) and Millennium 
(right). Each panel represents the midplane of the FP, as shown in Figure |2] The solid red line represents a least-squares fit to the 
binned data, p is the Spearman rank coefficient for the binned data, where p = 1.0 designates a perfect correlation, p = —1.0 designates 
a perfect anti-correlation, and p = 0.0 designates no correlation. Both SAMs demonstrate a strong correlation between age and velocity 
dispersion (p = 0.74 and 0.72 for SOS and Millennium, respectively), and a somewhat weaker anti-correlation between age and effective 
radius {p = —0.56 and -0.45 for SOS and Millennium, respectively). The Millennium metallicity has a very tight dependence on effective 
radius {p = 0.96) and a much weaker dependence on velocity dispersion (p = 0.39), while the SOS metallicity has a moderate dependence 
on both radius and velocity dispersion (p = O.Sl and 0.60, respectively). 



following a merger obeyed the same trends, with a smaller 
radius offset. As a result, the Millennium galaxies shown 
here have larger radii than observations would suggest. The 
SOS progenitors and remnants are, in general, more consis- 
tent with observations, though remnants with the lowest for- 
mation redshifts are systematically large by 0.2 dex. 

To determine the effect of the major merger itself on the 
age-FP and metallicity-FP relations, we have formed similar 
contours using the properties of the disk progenitors. Using 
the same methodology as described above, progenitors are 
placed along a plane relating surface brightness, effective ra- 



dius, and circular velocity. Burstein et al. (1997) found that 



the radii, surface brightnesses, and characteristic velocities 
of self- gravitating systems form a 'cosmic metaplane', where 
the slope of the respective plane varies with Hubble type. 
For our purposes, we we will call the late- type portion of 
this metaplane the late- type fundamental plane (LTFP). 

For consistency, we use the surface brightness that the 
progenitor galaxies would have had at redshift zero if they 
had not undergone the merger event and all star formation 
stopped at the time of the merger. While we do not expect 
this to be entirely accurate, we make this approximation 
in order to compare to the low-redshift results of G09. We 
note that LTFP and FP relations across a single slice are 
not directly related, since progenitors that lie below the FP 
may produce remnants that lie above it. The overall trends, 
however, are consistent across each LTFP slice, allowing us 
to make comparisons between trends in Re and Vdrc for the 
progenitors and those in Re and a for the remnants. 

Results are shown in Figures [9] and Overall, we find 
that progenitor ages increase with circular velocity and de- 



crease with radius, while their metallicities increase strongly 
with circular velocity and are nearly independent of radius. 

The differing age-FP and metallicity-FP relations for 
the progenitors stem from differences in the age-mass and 
metallicity-mass relations in the two models. Observations 
have shown that more massive spiral and elliptical galaxies 
tend to have older ages and higher metallicities than their 
less massive counterparts (Galla zzi et al.|2 005). While both 
the Millennium and SOS progenitors have increasing metal- 
licities as a function of mass, their ages are nearly indepen- 
dent of mass (Figure pT] see also Somerville et al. (2008a), 
Figure 6). Since the stellar mass is related to the circular 
velocity via the Tully-Fisher relation ( |Tully k Fisher|1977"| , 
M* oc V^j-c, a ^ 4, the age-mass and metallicity-mass cor- 
relations directly translate into age-LTFP and metallicity- 
LTFP relations for the progenitors; both age and metallicity 
should increase with circular velocity while remaining rela- 
tively independent of radius. 

The failure of the SAMS to reproduce the age-mass rela- 
tion can be seen clearly in Figure [9] as the radial dependence 
of the stellar age is in the opposite direction from what the 
age-mass relation would predict. Instead of a strong age- 
mass relation, we find that the progenitors' ages are linked 
to the redshift of the merger. Since we calculate all stellar 
properties at redshift zero but truncate star formation at 
the time of the merger, galaxies that merge at high redshifts 
have intrinsically older ages than those that merge at lower 
redshifts. This represents an artifact of the analysis meth- 
ods, which will be corrected when we include subsequent 
evolution in future work. As we will show in section 5, how- 
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Figure 6. Relation between central dark matter fraction, effective radius, and velocity dispersion for elliptical galaxies in SOS (top) and 
Millennium (bottom). The different panels represent slices of the FP, as shown in Figure [2] In both SAMs, galaxies that fall below the 
FP have higher central dark matter fractions. 



ever, the age-FP and metallicity-FP relations are similar 
regardless of the redshift of the merger. 

As the SAMs correctly reproduce the metallicity-mass 
relation for spiral galaxies, they predict a metallicity-LTFP 
relation that is strongly dependent on circular velocity and 



is nearly independent of radius (Figure 10). It is also worth 
noting that both the Millennium progenitors and remnants 
span a much smaller range in metallicity than their SOS 
counterparts. 

Comparing these relations with those of the remnants, 
it is apparent that there is a counterclockwise rotation in 
the Re — (J plane between the comparable relations for pro- 
genitors and remnants. This is made evident in Figure |12| 
in which we plot the age and metallicity FP and LTFP cor- 
relations from the midplane alongside each other. If we form 
a linear relation between the binned ages and metallicities, 
log Re, and \ogVchar, where Vchar = Vdrc for progenitors 
and (J for remnants, we find that the amount of rotation is 
:^ 20° counterclockwise for age and metallicity, respectively, 
in SOS. The amount of rotation in the age-FP correlation is 
slightly lower for Millennium (pi 10°), while the metallicity- 
FP correlation exhibits a higher amount of rotation (:^ 40°). 

This net rotation is due to the varying gas fractions of 
the progenitors. A galaxy's gas contents can be traced to its 
stellar mass, such that more massive galaxies tend to have 



lower gas fractions (Figure 10, see also Kannappan 2004 
Catinella et al.|["2010 ); furthermore, at fixed mass, the Ken- 



respondingly higher gas fractions (Covington et al. 2011). 
Both of these trends are evident when we use the progen- 
itors to produce gas fraction-LTFP and stellar mass-LTFP 



relations (Figure 13). The gas fraction is strongly correlated 
with the stellar mass, with the gas-surface density relation 
entering as a higher order effect. 

In the dissipational model, galaxies with higher gas frac- 
tions undergo more dissipation during the merger, producing 
remnants that are more compact and have lower velocity 



dispersions (see Section 2.2). Combining this with the de- 
pendence of gas fraction on stellar mass and surface density, 
we find that progenitors with low masses and low surface 
densities, which lie in the upper left of each of the LTFP 
slices, tend to contract the most, moving to smaller radii 
and larger velocity dispersions after the merger. Conversely, 
galaxies with high masses and high surface densities have the 
lowest gas fractions, producing remnants with larger radii 
and smaller velocity dispersions. Thus the rotation between 
the progenitor and remnant properties is largely an effect of 
gas fraction. Since the Millennium galaxies have lower gas 
fractions than the SOS galaxies, the amount of rotation Pre- 
vious studies have shown that a mass-dependent gas fraction 
can introduce the observed rotation between the virial plane 



and the fundamental plane (Dekel & Cox||2006 


Robertson 


et aL|2006 Hopkins et aL|200S| 


Covington et al. 


|2011); here 



nicutt relation implies that galaxies with larger radii have 
lower surface densities, less efficient star formation and cor- 



we point out the effects of that rotation on stellar population 
parameters. 

It is worth noting that the amount of star formation has 
a minimal impact on the ages and metallicities of the rem- 
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Figure 7. Relation between stellar mass-to-light ratio, effective radius, and velocity dispersion for elliptical galaxies in SOS (top) and 
Millennium (bottom). The different panels represent slices of the FP, as shown in Figure [2] Galaxies that lie above the FP tend to have 
lower mass-to-light ratios 
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Figure 8. Relation between stellar mass-to-light ratio and the redshift of the merger for galaxies in SOS (left) and Millennium (right). 
Different colors represent different slices of the FP, as shown in Figure |2] While galaxies that have very low merger redshifts tend to lie 
in the highest FP slice and have low mass-to-light ratios, there is no correlation between merger redshift and mass-to-light ratio in the 
other FP slices. 



nants; there is a close correspondence between the age and 
metallicity of the progenitor and remnant, with Spearman 
rank coefficients greater than p = 0.95 for both parame- 
ters and both SAMs. Thus, the merger has a larger impact 
on the structural properties of a galaxy than on its stellar 
population parameters. 



4.2 Comparison to observations 

To better compare with G09, we have 
plotted the age and metallicity contours over 
the range (1.9 km s~"^ < log a < 2.4 km s~^), 



re- 



(0.0 kpc < log Re < 0.7 kpc) considered by Graves et al 
(|2009b|) alongside the G09 data (Figures [m] and [15)" 



We caution that the G09 ages were later found to be 
systematically high by ^ 0.12 dex, owing to weak emission 
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Figure 9. Relation between luminosity-weighted age, effective radius, and circular velocity for spiral progenitor galaxies in SOS (top) and 
Millennium (bottom). The different panels represent slices of the late- type fundamental plane (LTFP). In both SAMs, stellar population 
age increases with circular velocity and decreases with radius. 



in the H/3 absorption line ( Graves &: Faber||2010| ; however, 
this would not effect the overall trends. G09 also calculated 
luminosity-weighted ages and metallicites using the Lick 
indices, which have been shown to more closely correlate 
with the period of most recent star formation, instead of a 
global quantity ( Trager Somerville^2009j . In future work, 
we will account for this discrepancy by measuring the Lick 
indices directly. 

Examining the trends within FP slices, the age-FP cor- 
relations are in rough agreement with G09. We caution that 
this agreement occurs in spite of the fact that the progen- 
itors fail to reproduce the observed mass-age relation. A 
stronger correlation between age and mass in the progeni- 
tors would directly lead to a stronger age- V^circ dependence 
in the LTFP, producing remnants with age-FP correlations 
that are dependent on both velocity dispersion and radius. 
Thus, improving the progenitors' fit to the mass-age relation 
of Gallazzi et aL' (2005) would have the effect of worsening 
the age-FP correlation when compared to G09. The major 
difference between our results and those of G09 is that we 
find metallicity to be dependent on radius and velocity dis- 
persion, while G09 found metallicity to be dependent on 
velocity dispersion alone. 

Both of the SAMs do a better job of reproducing ob- 
served trends through, as opposed to across, the FP. Galax- 
ies that fall above the FP tend to be younger and more 
metal-enhanced than average, while those that fall below 
the FP are older and more metal-poor, in agreement with 
G09. This is partially an artifact of truncating star forma- 



tion after the major merger: galaxies above the FP have 
universally low merger redshifts while those within and the 
below the FP have a wide range of merger redshifts. This re- 
sults in their having a more recent epoch of star formation 
and correspondingly younger stellar ages, as well as more 
extended star formation histories, leading to higher metal- 
licit ies. As discussed above however, even at fixed redshift, 
galaxies above the FP tend to have younger ages and higher 
metallicities, suggesting that the correlations through the 
FP are real and not merely artifacts of the methodology. 



5 EVOLUTION WITH REDSHIFT 

The assumption that elliptical remnants experience no fur- 
ther growth or other structural changes following their for- 
mation via a major merger has direct implications for the 
results. High-redshift observations point to a population of 
massive, compact ellipticals that are not seen in the local 



universe (Bell et al. 20M Trujillo et aL 2006^ Toft et al. 2007 



van Dokkum et al. 2008: Damjanov e t al.|2009| Franx et al.| 
2008| |van Dokkum et al. 2010; Wilh ams et al.|2010| ). There 
is growing evidence that subsequent minor dry mergers may 
account for at least some of this size evolution, particularly 
for the most massive galaxies ( Naab et al.||2009| Oser et al. 



2010 : Shankar et al.'2010"Newman et al.|2011||Shankar et al 



2011; Trujillo et al. 2011), but see also Nipoti et al. (2009); 



Nair et al. ( 2011 ); Saracco et al. ( 2011 ). Thus, we expect that 



the elliptical galaxies with high formation redshifts would 
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Figure 10. Relation between luminosity-weighted metallicity, effective radius, and circular velocity for spiral progenitor galaxies in SOS 
(top) and Millennium (bottom). The different panels represent slices of the LTFP. In both SAMs, metallicity increases with circular 
velocity but is nearly independent of radius. 
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Figure 11. From left: stellar gas fraction, age, and metallicity as a function of stellar mass for spiral progenitors from the SOS (blue 
diamonds) and Millennium (red diamonds) SAMs. The black circles represent SDSS observations from [Gallazzi et al.| ( [2005| >. In both 
SAMs, gas fraction and metallicity decrease with stellar mass in qualitative agreement with observations, though the Millennium pro- 
genitors span a much smaller range in metallicity than their SOS counterparts. However, both SAMs produce a flat age-mass relation, in 
disagreement with observations. 
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Figure 12. Relation between effective radius, characteristic velocity, age (left panels) and metallicity (right panels) for progenitor (top 
panels) and remnant (bottom panels) galaxies in SOS (top) and Millennium (bottom). All panels represent the middle slice of the FP 
and LTFP. In the progenitor panels, the characteristic velocity is Vdrc, while in the remnant panels the characteristic velocity is a. All 
of the relations show a counterclockwise rotation between the corresponding relations for the progenitors and remnants. The contours 
are individually normalized. 



undergo significant size evolution by z=0. The contours pre- 
sented so far stack all ellipticals in a slice, regardless of their 
formation timescale; here, we divide the population accord- 
ing to the redshift of the major merger when the ellipticals 
formed. This allows us to determine whether the redshift of 
the major merger has an effect on the FP correlations. Fig- 
ure [16] shows the age and metallicity FP correlations split 
into 3 redshift bins (0.0 < z < 1.0, 1.0 < z < 2.0, and 
2.0 < z < 3.0). 

Since galaxies that merge at low redshifts tend to be 
larger and more gas-poor, the low-redshift elliptical galaxies 
have both larger radii and smaller velocity dispersions on 
average. We find only minimal change in the age-FP trend, 
and no change in the metallicity-FP trend as a function of 
redshift. This is because the dominant factors in producing 
these trends, namely the relations between stellar mass and 
age, metallicity, and gas fraction, are similar at all redshifts 
(Figure 17). Galaxies that merge at lower redshifts tend to 
have younger ages and lower gas fractions, but the slopes of 
the mass-gas fraction and mass-age relations are relatively 
unchanged. 



5.1 Effects of subsequent evolution 

While the SOS and Millennium SAMs reproduce the ob- 
served elliptical mass function ( Somerville et al.|200Sb Cro- 
ton et a'L]|2006), the assumption that galaxies passively 



evolve following the major merger overproduces low-mass 
ellipticals and underproduces the high end of the mass func- 
tion. This is an expected artifact of the dissipational model: 
by prohibiting subsequent evolution, we preserve galaxies as 
they were when they first formed, preventing mass growth 
through star formation or mergers. Furthermore, we do not 



allow for bulge formation through disk instabilities. Overall, 
Cll found that applying the dissipational model to major 
mergers of disk-dominated galaxies alone reproduces only 
a third of the ellipticals seen in the local universe. Imple- 
menting the model directly within SAMs will reproduce the 
observed abundances, but this will require the model to be 
expanded to predict the velocity dispersion and radius after 
minor and mixed-morphology mergers. 

Although we have not attempted to model the evolution 
of ellipticals, we can estimate its effect. Dissipationless sim- 
ulations of major mergers of ellipticals suggest that, while 
the size-mass relation is preserved, there is increased scat- 
ter in radius due to variations in the orbits of the progen- 
itors (|Barn es"1992: Hernqui st et aLl|1993| |Cole et al.||2000' 



Shen et al.^2003 i Boylan-Kolchin et al.|2005t . This is gener- 
ally consistent with Cll, who found a size-mass dispersion 
that was lower than observations. Robertson et al. (2006); 



Boylan-Kolchin et al. ( 2006 ) demonstrated that major merg- 
ers of ellipticals do little to change the tilt of the FP; the 
increase in the effective radius is offset by similar adjust- 
ments in velocity dispersion and surface mass density. How- 
ever, recent SPH simulations of dry, minor mergers between 
a massive elliptical and disk galaxies have suggested that 
subsequent mergers will greatly increase the radius while 



slightly decreasing velocity dispersion (Naab et al. 2009 



Oser et ar]|2010| 2012). Thus, we suspect that subsequent 



mergers will weaken the correlations we have found between 
stellar population parameters and radius, without affecting 
the trends with velocity dispersion. This would bring our 
results into closer agreement with the analysis of G09, al- 
though the magnitude of the changes remains unclear. 

Alongside the binary major-merger scenario, additional 
pathways to elliptical galaxy formation have also been sug- 
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Figure 13. Relation between effective radius, circular velocity, gas fraction (top panels) and stellar mass (bottom panels) for spiral 
merger progenitor galaxies in SOS (top) and Millennium (bottom). The different panels represent slices of the LTFP. In both SAMs, the 
gas fraction is strongly correlated with the progenitor stellar mass; regions with high masses have low gas fractions, and vice versa. At 
fixed mass, galaxies with larger radii have lower surface densities and higher gas fractions. 



gested. "Naab et al." ('2007') have proposed minor mergers of 6 CONCLUSIONS 



disk galaxies, while Dekel et al. (2009) point to the impor- 



tance of "cold flows" of gas into disk galaxies. However, the 
current generation of SAMs is not equipped to simulate such 
events (Hirschmann et al.||2011|. We hope that including 
multiple minor mergers and cold flows in the future will im- 
prove the accuracy of our predictions. Both the Millennium 
and SOS SAMs also allow for bulge formation through disk 
instabilities, which we have also not attempted to model 
here. 



We have used initial conditions for disk-like merger progen- 
itors predicted by the SOS fSomerville et al. 200Sb) and 
Millennium ( |Croton et al.|[20 06) SAMs along with an an- 
alytic model for the size and velocity dispersion of merger 
remnants from Covington et al. (200S, 2011) to predict the 
distribution of these remnants across and through the fun- 
damental plane. We passively evolve the stellar populations 
and structural properties of these remnants to redshift zero, 
at which point we calculate their luminosity-weighted ages 
and metallicites and separate them according to their resid- 
ual surface brightness in the fundamental plane. 

For both SOS and Millennium, we find that elliptical 
ages increase as a strong function of velocity dispersion and 
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Figure 14. Relation between luminosity-weighted age, effective radius, and velocity dispersion for elliptical galaxies in G09 (top) SOS 
(middle) and Millennium (bottom). Here we plot only the region considered in G09. The different panels represent the three central 
slices of the FP, as shown in Figure [2] In the SAMs and the observations, stellar population age increases with velocity dispersion, but 
the simulated galaxies have a further dependence on radius. Galaxies that lie above the FP also tend to be younger than those that lie 
below the FP. 




Figure 15. Relation between luminosity-weighted metallicity, effective radius, and velocity dispersion for elliptical galaxies in SOS (top) 
and Millennium (bottom). Here we plot only the region considered in G09. The different panels represent the three central slices of the 
FP, as shown in Figure[2] While [Mg/H] depends strongly on velocity dispersion in G09, in the SAMs metallicity depends on both velocity 
dispersion and effective radius. The SAMs tend to have lower metallicities than observations on average; furthermore, the Millennium 
SAM produces galaxies with a very narrow range in metallicities. 



decrease with radius. The dependence on velocity dispersion 
is in agreement with G09 while the dependence on radius 
is not found in observations. The metallicities, however, in- 
crease with both radius and velocity dispersion, whereas the 
G09 metallicities are independent of radius. In all cases, the 
relations for ellipticals are rotated from the corresponding 
relations for the progenitors. This rotation arises because 



the amount of dissipation is strongly dependent on the gas 
fraction of the progenitor galaxies, which is correlated with 
the progenitors' masses and circular velocities. 

The difference between the age-FP and metallicity- 
FP correlations in the remnants stem from the different 
age-mass and metallicity-mass relations in the progenitors. 
While the metallicity-FP correlations are closer to obser- 



© 0000 RAS, MNRAS 000, 000-000 



18 Porter et al 



S08 Millennium 
0.0<z<1.0 1.0<z<2.0 2.0<z<3.0 0.0<z<1.0 1.0<z<2.0 2.0 < z < 3.0 




Log o [km/s] Log o [km/s] 

Figure 16. Relation between luminosity- weighted age (top) and metallicity (bottom), effective radius, and velocity dispersion for SOS 
(left) and Millennium (right), binned by the redshift of the major merger. Each panel represents the mid-FP, as shown in Figure [2] 
Horizontally, panels represent three formation epochs (0.0 < z < 1.0, 1.0 < z < 2.0, and 2.0 < z < 3.0); each galaxy is only plotted 
in the redshift bin corresponding to its major merger. The ages of each redshift panel are individually scaled; from left, they span log 
[0.61, 1.00], [0.95, 1.05],and[1.03, 1.07] Gyr for SOS and log [O.SO, 1.05], [0.99, 1.06], and[1.05, l.OS] Gyr for Millennium. Galaxies that merge 
at lower redshifts tend to be larger, with lower velocity dispersions. With the exception of the youngest galaxies within a redshift bin, 
neither the age-FP relation nor the metallicity-FP relation show an evolution with redshift. 




Figure 17. From left: Evolution of the progenitor stellar gas fraction-, age-, and metallicity-mass correlations for galaxies that merge 
between z = 3.0 and z = 0.0. SOS progenitors are shown in the top panels, while Millennium remnants are shown in the bottom panels. 
Although both SAMs show offsets such that galaxies that merge at higher redshifts have higher gas fractions and older ages at redshift 
0.0, the overall slopes of the relations do not vary with redshift. 



vations in the local Universe than the age-FP correlations, 
the spiral progenitors are uniformly too old, suggesting that 
the metallicity-FP correlation is closer to the 'true' correla- 
tion following a major merger. The differences between the 
results using simulated galaxies and the analysis of G09 sug- 
gest that a single major merger is not sufficient to remove the 
correlation between the progenitors' and remnants' radii. 

In agreement with G09, we find that galaxies with 
higher residual surface brightness tend to be younger and 
more metal-rich. Examining their structural properties, we 



find them to have lower stellar mass-to-light ratios and 
lower central dark matter fractions. These effects persist 
even when variations in mass-to-light ratio as a function 
of merger redshift are considered, suggesting that this is not 
merely a theoretical artifact. The galaxies with extremely 
high surface brightnesses do have lower formation redshifts, 
suggesting that they will fall to lower FP slices in the future, 
as their mass-to- light ratios increase. This is in accordance 



with [Forbes et al.| ( |1998| ; [Terlevich Forbes] ( |2002| ), and 
G09. 



© 0000 RAS, MNRAS 000, 000-000 



The Effect of Major Mergers on Age and Metallicity Across the Fundamental Plane 19 



We find that variations in the stellar mass- to- light ra- 
tio and the dark matter fraction within one effective radius 
contribute to the thickness of the FP, with the dark matter 
fraction having the larger effect. If we look across the FP, we 
find that the two quantities rotate the plane around different 
axes, with the mass- to-light ratio being highly dependent on 
velocity dispersion while the dark matter fraction is largely 
dependent on effective radius. 

The present framework is limited in that it assumes that 
elliptical galaxies are only formed via a single major merger 
between disk-dominated progenitors, and does not allow for 
any subsequent growth or mergers. As such, we only produce 
^ 1/3 of the ellipticals seen in the local universe, with the 
deficiency being strongest at the high-mass end of the mass 
function. In future work, we will incorporate a broader range 
of processes, including subsequent minor mergers, elliptical- 
elliptical mergers, and mixed morphology mergers. We will 
also separate the red sequence population considered in G09 
into elliptical and SO/Sa galaxies, using an automated mor- 
phological classification scheme ( [Cheng et al.||2011| . 

A new version of the Somerville SAM also contains a 
detailed Galactic Chemical Evolution model, including non- 
instantaneous recycling, enrichment by both core collapse 
and Type la supernovae, and tracking of multiple chemical 

( |2010b|a| ). We also 



elements, as described in Arrigoni et al, 



plan to calculate ages and metallicities using the Lick indices 
instead of using approximations to calculate the luminosity- 
weighted properties. This will allow us to make direct com- 
parisons with [Fe/H], [Mg/H], and [Fe/Mg]. These improve- 
ments will allow us to make more detailed comparisons to 
observations. 
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